function y = func_n1_from_p(p, CONSTS)
    eta = CONSTS.eta;
    n2 = func_n2_from_p(p, CONSTS);
    
    y = - (eta./n2);
end
